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We present a theoretical study of the influence of resonant enhancement on quantum path dynam¬ 
ics in the generation of harmonics above and below the ionization threshold in helium. By varying 
the wavelength and intensity of the driving held from 425 nm to 500 nm and from 30 TW/cm 2 
to 140 TW/cm 2 , respectively, we identify enhancements of harmonics 7, 9, and 11 that correspond 
to multiphoton resonances between the ground state and the Stark shifted ls2p, ls3p, and ls4p 
excited states. A time-frequency analysis of the emission shows that both the short and long quan¬ 
tum path contributions to the harmonic yield are enhanced through these bound state resonances. 

We analyze the sub-cycle time structure of the 9th harmonic yield in the vicinity of the resonances 
and find that on resonance the long trajectory contribution is phase shifted by approximately 7t/4. 

Finally, we compare the single atom and the macroscopic response of a helium gas and find that 
while the sub-cycle time profiles are slightly distorted by propagation effects, the phase shift of the 
long trajectory contribution is still recognizable. 

PACS numbers: 42.50.Hz, 42.65.Ky, 32.80.Rm 


I. INTRODUCTION 

High harmonic generation is a versatile source of ultra¬ 
fast, coherent extreme ultraviolet (XUV) radiation i.l/2j, 
produced by the interaction between an intense infrared 
or visible laser held and a gas of atoms or molecules. 
The semi-classical understanding of harmonic generation, 
in which an electron wave packet is initiated by tunnel 
ionization, accelerates in the laser held, and finally re¬ 
scatters on the parent ion mm has given rise to the 
held of high harmonic spectroscopy (HHS), in which the 
rescattering electron wave packet has been used as a sen¬ 
sitive probe of the structure and dynamics of the ion core 

mm- 

Resonantly enhanced high harmonic generation (RE- 
HHG) is a particular example of the more general HHS 
technique. Although REHHG has been shown to increase 
the harmonic yield in a limited range of experimental 
settings HDHUl, it has more recently been explored for 
its potential to learn about the dynamics of bound and 
quasi-bound states in the presence of a strong driving 
held [J3J ITO|2;)I . Several mechanisms for resonant en¬ 
hancement have been discussed in the literature, gener¬ 
ally all involving an intermediate, resonant, step in the 
semi-classical model. The resonant step may occur ei¬ 
ther in the ionization process, via a multiphoton reso¬ 
nance between the ground state and the Stark-shifted 
excited state, or in the rescattering process via enhanced 
recombination, or by capture into an excited bound state 
that subsequently decays via spontaneous emission of 
light P2, I! [ TH OS [20] dT • The capture and sponta¬ 
neous emission process has been explored in detail for 
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short-lived quasi-bound states embedded in a continuum 
for which it can give rise to very large enhancements 
Hang. For bound-state resonances with long lifetimes, 
the capture and spontaneous emission process can gener¬ 
ally be distinguished from the coherently driven resonant 
enhancement (via multiphoton ionization or enhanced re¬ 
combination), since it will give rise to narrow-band ra¬ 
diation at the field-free resonance frequency given that 
it largely takes place after the driving laser pulse is over 
33 3H- br contrast to this, the coherently driven res¬ 
onantly enhanced response will give rise to emission at 
the difference frequency between the ground state and the 
Stark-shifted excited state since this process only takes 
place while the laser field is on. 

In this paper we study the coherently driven REHHG 
and investigate the interplay between the resonant en¬ 
hancement and the quantum path dynamics of the har¬ 
monic generation process. In particular, we study how 
the amplitude and phase of the different quantum path 
contributions to the harmonic yield in helium are changed 
in the vicinity of a bound state resonance. We do this 
by first solving the time-dependent Schrodinger equation 
(TDSE) in the single active electron (SAE) approxima¬ 
tion for wavelengths between 425 — 500 nm and for inten¬ 
sities up to 140 TW/cm 2 . We have chosen the intensity 
and wavelength regime such that the Keldysh parame¬ 
ter, of approximately 2, is in line with previous studies 
PUQSIEI]. We identify resonant enhancement of har¬ 
monics 7, 9, and 11 (H7, H9, Hll) due to Stark-shifted 
resonances between the ground state and the ls2p, ls3p, 
and IsAp states. By considering the time-frequency pro¬ 
files of the harmonic emission, we also clearly identify 
the familiar short and long quantum path contributions 
to these harmonics both above and below the field-free 
ionization threshold. We find that close to the ls3p reso¬ 
nance, the amplitude of both the short and long quantum 
path contributions are enhanced. We also find that the 
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long-trajectory contribution is phase shifted by approx¬ 
imately 7t/4 in the vicinity of this resonance, while the 
short-trajectory contribution is not. These results indi¬ 
cate that while the enhancement of the yield happens in 
the ionization step, the phase shift is due to the inter¬ 
action between the returning electron wave packet and 
the atomic potential, for which there is a large difference 
between the short and long trajectory dynamics. Finally, 
we calculate the macroscopic harmonic response by solv¬ 
ing the coupled Maxwell wave equation and the TDSE 
for a gas of atoms interacting with a focused laser beam. 
We find that the enhancement and the phase shift of the 
resonantly enhanced harmonics can still be recognized in 
the macroscopic response. 

The paper is structured as follows: Section II presents 
a brief introduction to the theoretical methods used in 
this paper. Section III and IV focus on the single atom 
response; in Section III we look at REHHG in the fre¬ 
quency domain and use the harmonic spectra to map out 
the enhancement as a function of driving intensity and 
wavelength, and in Section IV we examine REHHG in 
the time domain and study the effect of REHHG on the 
quantum trajectories. Section V presents results of our 
macroscopic calculations, and Section IV presents a sum¬ 
mary of our findings. 

II. THEORETICAL METHOD 

We numerically solve the TDSE in the SAE approxi¬ 
mation for a He atom interacting with an intense laser 
field. We use a pseudo-potential that reproduces well 
the energies and oscillator strengths of the singly excited 
states of helium to describe the interaction of the active 
electron with the held. We will generally refer to the 
excited 1 snp states as np. For more details on our nu¬ 
merical method and the pseudo-potential see [25]. The 
laser held has the form of a cosine squared pulse with a 
full width at half maximum (FWHM) duration r, given 
by: 

E L (t) = E 0 cos( —) 2 cos(27rf), (1) 

T 

where Eq is the peak electric held, C 2 = 2 arccos(| 4 ), and 
t is time in units of optical cycles (O.C.) of the driving 
held [55]. With a FWHM duration r = 6 O.C. the total 
duration of the pulse, and thereby the TDSE integration, 
is 16.5 O.C. We calculate the dipole spectrum from the 
time-dependent acceleration a(t) as D(cj) = —A(oj)/uj 2 : 
where the acceleration spectrum is given by: 

A(u) = J a{t)W H {t)e iut dt. (2) 

In this equation, Wn{t) = -5(1 + cos(^)) is a Hanning 
window function |27j which brings the time-dependent 
acceleration smoothly to zero at the end of the calcu¬ 
lation. In numerical calculations it is always necessary 


to include a window function on a(t ) when there is any 
population remaining in excited states at the end of the 
pulse, since the coherence between the excited and the 
ground state population gives rise to a dipole moment 
which in the calculation would last forever. In an ex¬ 
periment, this would give rise to spontaneous emission, 
decaying with the life time of the excited states. The win¬ 
dow function thus serves to impose an artificial lifetime 
on the excited states. In most of the calculations shown 
in this paper, we use a Hanning window function with 
th = 8.25 O.C. in order to last only as long as the driv¬ 
ing laser pulse. This means we are primarily considering 
the coherent, driven, response of the atom to the laser 
pulse and not the radiation produced by long-lived decay 
of excited state populations. It also means that in most 
of the calculations shown, we are discarding the capture 
and spontaneous emission contribution to the resonantly 
enhanced response. By doing a few calculations where we 
continue the TDSE integration after the end of the laser 
pulse (using a longer Hanning window where th = 8.25 
O.C. for t < 0 and th = 30 O.C. for t > 0) we will show 
that this latter contribution gives rise to narrow peaks at 
the position of the field-free resonant frequencies which, 
at the intensities used in this study, is generally outside 
or in the wings of the harmonic spectral bandwidths. 

We also calculate the macroscopic harmonic response 
by solving the coupled Maxwell wave equation and TDSE 
(MWE-TDSE) via space-marching of the full electric 
field E(oj), containing both laser and harmonics fields, 
through the nonlinear medium [281 . In the slowly evolv¬ 
ing wave approximation, we can express the MWE in the 
following form (in SI units): 

C) • r\ 2 

Vi E{u) + —^-E{uj) = —^[P(w) + hon M], (3) 
c oz eoc 

where E(w), P(w), and Pj on (w) are also functions of the 
cylindrical coordinates r and z. P(w) = 2 N at D(uj) is the 
macroscopic polarization held, where N at is the atomic 
density and D(u>) is the one-electron single-atom dipole 
moment. Pi on is related to the free-electron contribution 
to the refractive index and is ignored for these calcu¬ 
lations as it is a small term for the parameters we are 
studying. At each step in the propagation direction z, 
we calculate P(w) in Eq. ^ by solving the TDSE at 
each radial point, and then propagate Eq. ([3| one step to 
get to the next z plane. More details of our MWE-TDSE 
solution can be found in [55J152]. 

In the macroscopic calculations, the focused laser beam 
is modeled as a Gaussian beam with a peak intensity of 
140 TW/cm 2 , a confocal parameter of 1 cm, and the 
focus of the beam is 1 mm before the center of the 1 mm 
long gas medium that has a pressure of approximately 30 
Torr. We use the same temporal profile as in the single 
atom calculations. All macroscopic calculations shown in 
this paper use the radially integrated yield of the electric 
field intensity at the end of the medium. 
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FIG. 1. Harmonic spectra for different driving wavelengths 
for a driving laser intensity of 140 TW/cm 2 . 



III. SINGLE ATOM: FREQUENCY DOMAIN. 

In this section we characterize the resonant enhance¬ 
ments in the harmonic spectrum as a function of laser 
wavelength and intensity. By ’’resonant enhancement” 
we mean that, for a given intensity, a harmonic is much 
stronger at a particular wavelength than at others be¬ 
cause the harmonic energy is resonant with a transi¬ 
tion between the ground state and a Stark-shifted ex¬ 
cited state GSBEH. For most of the np-states, the Stark 
shift is approximately equal to the ponderomotive shift of 
the continuum states mi mi mi- Thus direct resonant 
enhancement of harmonic q occurs when the following 
equation is satisfied : 

I E np — Eq \ + Up = qhui , (4) 

where E np and Eq are the field-free energies of the up¬ 
state and the ground state of helium, respectively, U p 
is the ponderomotive energy, and Hco is the laser pho¬ 
ton energy. Indirect resonant enhancements occur when 
Eq. (4) is satisfied and a harmonic q' near harmonic q is 
enhanced. 

Fig. [I] shows harmonic spectra for driving wavelengths 
of 426, 440, and 483 nm at an intensity of 140 TW/cm 2 . 
We see two enhancements in the spectra, H9 at 440 nm 
and H7 at 426 nm. Eq. (4) shows that for these three 
wavelengths, it is only the enhancement of H9 at 440 nm 
which represents a direct enhancement, via a 9-photon 
resonance between the ground state and the Stark-shifted 
3 p state. As we will show below, the enhancement of H7 
at 426 nm is an indirect enhancement, via a 9-photon 
resonance between the ground state and the 4 p state. 

Fig. [2] shows the spectral yield for a range of driving 
wavelengths between 420 nm and 500 nm. We have used 
two different Hanning windows to be able to distinguish 
between the coherently driven response, and the response 
due to the capture and spontaneous emission process. 


FIG. 2. Comparison of the harmonic spectra as functions of 
driving wavelength and energy for different Hanning windows 
for a driving intensity of 140 TW/cm 2 . (a) shows the spectra 
for a Hanning window that matches the duration of the pulse, 
(b) shows the spectra for a Hanning window that allows the 
dipole to ring for approximately 20 optical cycles after the 
pulse ends. The white dotted lines mark the position of H7- 
Hll as a function of driving wavelength. In (b) we have also 
indicated the field-free energies of the 2p-5p states. We note 
that the harmonic yield is shown on a linear scale, in arbitrary 
units. 


The shorter window in Fig. |2ja) shows the coherently 
driven response that we are predominantly interested in, 
whereas the longer window used in Fig. §b) addition¬ 
ally shows the radiation due to population left in excited 
states at the end of the pulse. The narrow spectral lines 
that are visible at the field-free transition frequencies in 
Fig. [2^b) are useful in identifying which state is primar¬ 
ily responsible for a given resonant enhancement. For 
example, the enhancement of H9 at 440 nm (470 nm) 
is clearly associated with a large population transfer to 
the 3 p (2 p) state whereas the modest enhancement of H7 
at 426 nm is associated with 9-photon-driven population 
transfer to the 4 p state. Interestingly, this 4 p resonance 
does not lead to a particular enhancement of H9, at least 
in comparison with the large enhancement due to the 3 p 
state. It is worth noting the large spectral separation 
between the coherently enhanced radiation and the field- 
free peaks due to the long-lived dipole moment. 

The resonant enhancements in any given harmonic can 
also be followed by considering the combined wavelength 
and intensity dependence of that harmonic. We start by 
showing in Fig. [3] the wavelength dependence of the H9 
yield, using an intensity of 140 TW/cm 2 . This plot cor¬ 
responds to a line-out of Fig. [2](a) along the white dashed 
line indicating the central frequency of H9, spectrally 
integrated from H8 to H10. At this intensity, one can 
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FIG. 3. H9 spectrally integrated yield as a function of driving 
wavelength, for an intensity of 140 TW/cm 2 . We show both 
the single atom yield (solid diamonds) and the macroscopic 
yield (open triangles). We have used the shorter Hanning 
window in calculating the harmonic spectrum. 


recognize two enhancement features in Fig. [3j a strong 
enhancement centered at 440 nm and a weaker enhance¬ 
ment centered at 476 nm, as we would expect from Fig. 
[2] The location of these peaks agree well with the pre¬ 
diction of Eq. (4) for the 3 p and 2 p states and q = 9. 

Next, Fig. [4] shows the combined intensity and wave¬ 
length dependence of H9, by plotting the H9 wavelength 
dependence as shown in Fig. [3] for many different driving 
intensities. For each intensity, the wavelength dependent 
yield has been normalized to allow for direct comparison 
between high and low intensities, whose true yields differ 
by orders of magnitude. The white lines in the figure 
indicate the predicted H9 photon energies where direct 
resonant enhancement would occur according to Eq. (4). 
Two enhancement features can be recognized in Fig. [4j 
due to the 3 p and 2 p states. The 2 p enhancement in gen¬ 
eral follows Eq. (4), but can be seen to shift less than 
ponderomotively for higher intensities. The 3 p enhance¬ 
ment feature follows the prediction of Eq. (4) only for 
intensities above approximately 60 TW/cm 2 . At lower 
intensities, the 3 p feature splits into two, with the lower 
energy branch marked by the dashed black line. 

The split 3 p enhancement feature can be understood 
as a generalized Autler-Townes splitting due to a near¬ 
resonant coupling between the 3 p and 2s states induced 
by the driving laser field. In this case, the dressed states 
of the laser-driven atom constitutes an Autler-Townes 
doublet which is approximately symmetric around the 
3 p state (32]. We note that the separation of the two 
’’states” appearing in the 3 p enhancement feature is not 
exactly equal to the Autler-Townes energy separation one 
would see in a fixed-pump, scanning-probe scenario. This 
is because in our calculation the pump (the driving laser) 
and the probe (the 9th harmonic) are locked to each other 
and thus both vary, which means that each of the dressed 
states in the enhancement map is 9-photon resonant with 
the ground state at a different pump wavelength. We be- 
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FIG. 4. H9 spectrally integrated yield as a function of driving 
wavelength and intensity. For each intensity, the spectrum 
has been normalized. The white lines represent the photon 
energies of the resonant enhancements due to the 2 p and 3p 
states predicted by Eq. (4). The dark dotted line highlights 
the lower energy feature of an Autler-Townes splitting of the 
3p feature (see text). 

lieve this is the first time an Autler-Townes splitting has 
been characterized using only the harmonic spectrum. 

To summarize this section: we have examined resonant 
enhancements in the harmonic spectrum using several 
different representations. We have shown that we can 
consistently identify and follow resonant enhancements 
in several harmonics due to the (Stark-shifted) 2 p — 4p 
states. We also showed that the 3 p enhancement feature 
splits into an Autler-Townes like doublet when the driv¬ 
ing laser wavelength is such that the 3 p is near-resonant 
with the 2s state. 


IV. SINGLE ATOM: TIME DOMAIN 

In this section we concentrate on the time profile of the 
harmonic emission in the vicinity of resonant enhance¬ 
ment. We will study both the overall envelope of the 
harmonic pulse, and its sub-cycle time structure. We cal¬ 
culate the time profile by selecting a range of frequencies 
from the harmonic spectrum and inverse Fourier trans¬ 
forming to the time domain (331 GH] ■ We start by show¬ 
ing the time-profile of the envelope of H9, calculated by 
applying a narrow (2w wide) frequency filter. Fig. [ 5 ] 
compares the time profile of 3p-enhanced H9 at 440 nm 
and 140 TW/cm 2 using the short and long Hanning win¬ 
dows, respectively. This comparison shows that the short 
Hanning window makes very little difference to the driven 
part of the harmonic response, but that it effectively cuts 
off the population-driven dipole response, in agreement 
with what we found in the spectral domain in Fig. [2] 
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FIG. 5. H9 time profile at 440 nm and 140 TW/cm 2 with the 
(scaled) intensity envelope of the driving pulse, for the short 
and long Hanning windows. 

We note that the oscillatory behavior in the tail of the 
H9 time profile is due to quantum beating between pop¬ 
ulation left in the 3 p and Ap states which due to their 
energy difference of 0.7 eV gives rise to a beat period of 
6 fs. Quantum beating as a result of a broadband XUV 
excitation has recently been experimentally observed in 
neon in a transient absorption scenario [35j . 

We next consider the sub-cycle time profile of the har¬ 
monic radiation by using a wide (12a;) spectral window. 
By sliding the central frequency of the window function 
through the harmonic spectrum, we construct the time- 
frequency profile of the harmonic radiation. An example 
of such a time-frequency profile is shown in Fig. §a) for 
the on-resonance case of driving wavelength and inten¬ 
sity of 440 nm and 140 TW/cm 2 . The sub-cycle time- 
frequency profile clearly shows that the harmonic gen¬ 
eration process, even for these low-order harmonics, is 
dominated by the well-known three-step model dynam¬ 
ics which leads to two dominant emission peaks per half¬ 
cycle, with the peaks getting closer in time for higher 
order harmonics. This is in agreement with previous 
findings pS I36H39] . Note that although the field-free 
ionization threshold for this wavelength is just below H9, 
the Stark-shifted ionization threshold at the peak of the 
pulse is well above H9, and the 9tli harmonic is thus a 
below-threshold harmonic for most of the duration of the 
pulse. 

Fig. |6](b) shows line-outs centered on H9 in the time- 
frequency plot, for three different driving wavelengths 
(scaled for better comparison). The long and short tra¬ 
jectory peaks (marked in the figure) have been identi¬ 
fied from time-frequency profiles such as that shown in 
Fig.§ a). For the short trajectory contribution, the fig¬ 
ure shows that as the wavelength increases the emission 
peak shifts to earlier in the cycle. This is as one would 
expect from the positive atto-chirp of the short trajec¬ 
tory contribution (40]. The negative atto-chirp of the 
long-trajectory contribution predicts that the peak would 
shift later in the cycle as the wavelength increases, which 



FIG. 6. (a) Time-frequency profile of the harmonic radia¬ 
tion for a driving wavelength of 440 nm and intensity of 140 
TW/cm 2 , when H9 is resonant with the Stark-shifted 3 p state, 
(b) Line-outs of the time-profile of H9 for driving wavelengths 
of 426 nm, 440 nm, and 483 nm. Note that the long and 
short trajectory emission peaks in (b) belong to successive 
half-cycles of the driving field as can be seen in (a). 


is true for the two non-resonant wavelengths. However, 
on resonance the long trajectory peak is shifted sharply 
later in the cycle by approximately 1/8 O.C, correspond¬ 
ing to a phase shift of 7r/4. 

We can follow the phase shift of the long trajectory 
peak near resonance by plotting the H9 sub-cycle time 
profile as a function of driving wavelength, as shown in 
Fig. □ The short and long trajectory peaks have again 
been identified from time-frequency profiles at different 
driving wavelengths. The white lines have been added 
as an approximate measure of the atto-chirp, by con¬ 
necting the emission times at the shortest and longest 
wavelengths that are both non-resonant (note that the 
traditional atto-chirp is not defined for the harmonics 
below the held-free threshold). As the wavelength in¬ 
creases from approximately 425 nm to 440 nm, the long 
trajectory emission peak shifts later in time by approx¬ 
imately 1/8 O.C. Between the 3 p and 2 p resonances at 
440 nm and 475 nm, respectively, the harmonic yield is 
weak and there are more than two emissions peaks per 
half O.C., which makes it difficult to identify a short and 
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FIG. 7. Sub-cycle time profile of the radiation centered on 
H9 as a function of driving wavelengths spanning the 3p and 
2p resonant enhancement features, which are located at 440 
nm and 470 nm for this intensity. The white lines are an 
approximate measure of the expected emission times for the 
short and long trajectory emission peaks. 


long trajectory contribution. However, as the wavelength 
increases toward the 2 p resonance at 475 nm, the short 
and long trajectory emission peaks are again clearly iden¬ 
tifiable. Although we cannot follow the quantum paths 
in the wavelength region below the 2 p resonance at 475 
nm, the shift of the long trajectory peak appears to be 
smaller than for the 3 p resonance. 

We find similar results for Hll at longer wavelengths 
(not shown in the figure), when it is in resonance with 
the 2p—5p states. As the (Stark-shifted) 4 p and 5 p states 
come into resonance around 490 nm, the long trajectory 
emission peak is shifted later in time by approximately 
1/10 O.C and stays shifted through the 3 p resonant en¬ 
hancement at 505 nm. For wavelengths between the 3 p 
and 2 p resonances (515-530 nm), the time profile is domi¬ 
nated by emission peaks that are not easily identifiable as 
long or short trajectory contributions, but as Hll comes 
into resonance with the Stark-shifted 2 p state around 540 
nm, these peaks reappear clearly in the time profile. 

The results in Figs. [6] and [T] allow us to draw a num¬ 
ber of conclusions: (i) Harmonics that are enhanced by 
bound state resonances still exhibit strong features that 
correlate with the well-known short and long re-colliding 
electron trajectories. Our results in fact suggest that 
the short and long trajectory contributions are more eas¬ 
ily identified close to resonance than in between, (ii) 
Both quantum path contributions are enhanced around 
the resonances, indicating that the enhancement happens 
predominantly in the ionization stage, via a multiphoton 
resonance between the driving field and a Stark-shifted 
resonance, (iii) The fact that only the long trajectory 
emission peak is shifted in time suggests that this phase 
shift is imposed in the recollision stage of the harmonic 


generation process - whereas the short trajectory emis¬ 
sion happens very shortly after the ionization time, the 
long trajectory wave packet spends a long time in the 
continuum and samples the full spatial range of the ionic 
potential upon return. Our results also suggest that this 
phase shift is specific to a particular resonance (the 3 p 
state), since both H9 and Hll show exhibit similar phase 
shifts around the 3 p resonance and almost no phase shift 
around the 2 p resonance. 

To conclude this section on the single-atom emission 
times, it is interesting to note that on resonance, the en¬ 
velope of the H9 emission is shifted slightly later in time 
than off resonance. This can be seen both in Fig. |6ja) , 
where the most intense pair of short/long emission peaks 
is centered around 1.2 O.C, and in Fig.[5]as the position 
of the envelope of the H9-only radiation which also peaks 
around 1.2 O.C. Off resonance, the envelope peaks closer 
to 0.7 O.C. (not shown in the figure). We also see this 
shift of the envelope for Hll around the np resonances. 
It is possible that this overall delay could be due to trap¬ 
ping of the electron in the excited state prior to ionization 
(delayed ionization), so that the excited electrons are still 
launched into the continuum at the times that lead to re¬ 
turn along the short or long trajectories, respectively, but 
half a cycle later. Such a scenario was discussed in [44] in 
the context of double ionization, in which an excited elec¬ 
tron is caught in a so-called sticky region characterized 
by resonances in the combined laser-ion potential. We 
note that it is not likely that the overall delay is caused 
by the electron being trapped in the excited state upon 
returning to the vicinity of the core - this would generally 
give rise to emission times that would not correspond to 
those of the short and long trajectory return times. 


V. MACROSCOPIC RESPONSE 

Finally, we consider the macroscopic harmonic re¬ 
sponse, in particular whether and how the resonant en¬ 
hancement manifests itself in this respect. The open tri¬ 
angles in Fig. [3] show the macroscopic yield of H9 as 
a function of wavelength, in direct comparison with the 
single-atom result (solid diamonds). The shape of the 
macroscopic curve agrees well with the shape of the sin¬ 
gle atom curve, the macroscopic results reproducing the 
strong 3 p enhancement feature and the weaker 2 p en¬ 
hancement feature, however, the macroscopic enhance¬ 
ment features are shifted slightly to longer wavelengths. 
This difference between the single atom and macroscopic 
calculations is a consequence of the focal volume aver¬ 
aging of the intensity that is present in the macroscopic 
calculation. Looking at Fig. [4j we see that at lower 
driving intensity, the resonant wavelength for both the 
3 p and 2 p single-atom enhancement features shifts to¬ 
wards longer wavelengths due to the change in the Stark 
shift. Overall, we conclude from Fig. [3] that the enhance¬ 
ments found at the single atom level do indeed survive 
the macroscopic propagation. 












7 



FIG. 8. Sub-cycle time profiles for the macroscopic H9 yield, 
for three different driving wavelengths. 

In Fig. [8j we compare the macroscopic sub-cycle time 
profiles centered on H9 for three different driving wave¬ 
lengths. We have slightly shifted the driving wavelengths 
from the single atom results of Fig. ®b) from 440 nm 
to 444 nm and 483 nm to 490 nm in order to compen¬ 
sate for the shift of the enhancement peaks to longer 
wavelengths in our macroscopic calculation as discussed 
above. We find good agreement between Fig. [8] and Fig. 
®b). While there are some distortions in the macro¬ 
scopic case, the long and short trajectory peaks for all 
driving wavelengths of Fig. [6]Jb) are recognizable in Fig. 
[8] at the corresponding shifted wavelength. The phase 
shift of the long trajectory peak is also still clearly rec¬ 
ognizable in the macroscopic response. One difference 
between the single atom and the macroscopic results is 
the relative enhancement of the long and short trajectory 
contributions. In the single atom case, both quantum 
path contributions are equally enhanced, whereas in the 
macroscopic case, the long trajectory contribution is rel¬ 
atively stronger, by approximately a factor of two. This 
can be understood by looking at the single atom calcu¬ 
lations in Fig. [7] For intensities slightly lower than the 
resonant peak intensity, the long trajectory contribution 
is already both enhanced and shifted compared to the 
short trajectory contribution. Since the macroscopic sig¬ 
nal results from the radially integrated harmonic yield 
(i.e. an integration over lower-intensity contributions), 
this favors the long trajectory contribution. Note that 
we have tested that this effect is not dependent on phase 
matching. The relative position of the laser focus and the 
center of the gas jet was chosen so as to optimize phase 
matching of the short trajectory, but in this case we in 
fact find little difference in the macroscopic yield or time 
profiles when we move the focus to the center of the jet. 


VI. SUMMARY 

We have presented a study of the interplay between 
resonant enhancement and quantum path dynamics in 
near-threshold harmonic generation in helium. We con¬ 
centrated on the driven harmonic-generation response by 
time-filtering the harmonic signal so as to suppress the 
long-lasting radiation that would result from population 
left in excited states at the end of the pulse. By varying 
the wavelength and intensity of the near-visible driving 
laser field, we have identified a number of direct and in¬ 
direct enhancements of H7, H9, and Hll via the Stark- 
shifted 2 p — 5 p states. For H9, we observed an Autler- 
Townes-like splitting of the enhancement feature due to 
the 3 p state, when the wavelength and intensity are such 
that the driving field strongly couples the 3 p state to the 
nearby dark 2s state. 

In terms of the quantum path dynamics, we found that 
both the short and long trajectory contributions to the 
harmonic emission can be easily identified for harmon¬ 
ics that are resonantly enhanced via the Stark-shifted np 
states. We found that both contributions are enhanced 
on resonance, and that the maximum of the envelope 
of the resonant harmonic is delayed by approximately 
0.5 O.C. We interpret this to mean that the enhance¬ 
ment happens via a multiphoton resonance between the 
ground state and the Stark-shifted excited state, and that 
the electron is then trapped for a while in the excited 
state before entering the continuum. Furthermore, we 
found that only the long trajectory contribution acquires 
a phase shift, which leads to a delay in emission time 
of approximately 0.125 O.C, suggesting that the phase 
shift is acquired in the interaction between the returning 
electron wave packet and the ion core for which there 
is a large difference in the short and long trajectory dy¬ 
namics. Finally, we showed that both the enhancement 
and the phase shift are still visible in the macroscopic 
response. This means that these effects could potentially 
be explored experimentally, especially considering that 
our calculations predict that the macroscopic response 
is dominated by the long trajectory contribution which 
exhibits the on-resonance phase shift. 
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